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SUMMARY 


An improved computer program was devised to calculate the average 
blast impulse loads acting on a wail of a cubicle when an explosive 
charge is detonated within the cubicle. It was formulated by the AED's 
Process Fngineering Laboratory in connection with the Safety Design 
Criteria Program. 


PROGRAM DFSCRIPTION 


The original computer proyram was to calcuiate data points for 
the construction of impulse charts in regulatory Department of Defense 
publication Technical Manual 5-1300, "Structures to Prevent the Effects 
of Accidental Expiosion." This original program was prepar?d by Acmann 
& Whitney, a consulting engineering firm in New York City, dated June 1969, 
under contract to Picatinny Arsenal. 


The input to the program was simplified. Instead of a possible 
five data cards-per—problem only two data cards were required. Two 
subroutines were added to calculate geometrical ratios and to specify 
the grid system needed in the problem solution. The output was modified 
to echo the input and give a clear print-out of the reflected impulse 
on each wall and the totai impulse acting on the wall in question. This 
program is suitable for use by an engineer with little or no computer 
background. 


The program can be used as a supplement to ™ -5-1300. It eliminates 
the necessity of interpolation or extrapolation from the impulse charts 
in the maouval and can save many hours of tedious hand calculations by 
completing wulciple impuzse caiculations in a few minutes. 


Step-by-step procedures are given ter specifying the computer input. 
Sample input sheets, problem solutions by computer and impulse charts of 
™ 5-1300 arc in Appendix A. 


The program prints the solution of each problem on 3 separate page. 
The output consists of the citle, input data, and the calculated impulse 
for cach reflecting surface and total impulse load. 

The computer program is written in Fortran IV and has been run on 


an IBM 360, Model 65. A copy of the Fortran coding and required input 
data is in Appendix B. 


Limitation 


Because of the limitations of the test data used to develop this 


computer program, extrapolation beyond this range may give inaccurate results. 


However, to overcome this, the restrictions of the geometrical ratios 
mentioned on Pages 4-12 of TM S-1300 are incorporated in the computed 
progran. 
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INPUT IEFINITIONS AND FIELD SPECIFICATIONS 


Card No. 1 


Card No. 2 


Card No. 2. 11-20 W 
Card No. 2 31-1,0 ‘ e 
Card No. 2 lL-50 f 4 
cA 
Card Xo. 2 61-60/" 4 
ff 
ie 


Card No. 2 61-70 


Card No. 2 71 
72 
73 
7h 
\ 
eo 
ra 


TABLE ] ; 


f 


ExploZive charge weight (1bs.) 
(Decdmal point required.) 


/ 
sight of the wall in question (ft.) 
(Decimal point required.) 


Length of loaded tall (ft.) 
(Decimal point required.) 


Height of charge above floor (ft.) 
(Decimal point required.) 


Minimum distance between the charge 
to and an adjacent. wall (ft.) 
(Decimal point required.) 


Leave blank 


FLOOR - Insert 1 if prssent 
ROOF - Insert 1 if  Uresent 
LEFT SIDE - Insert 1 i: present 
RIGHT SIDE- Insert 1 if present 
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STEP-BY-STEP PROCEDURE FOR PROGRAM INPUT 


To utilize the program, it is necessary to fill out the information 
required by the two data cards on the input form (Figure 4). The data 
should be key-punched and inserted in its proper place in the input deck. 
It should be noted that the required input data in Appendix B must 
precede this data. 


The input parameters are given by the input sheet (Figure 4) and 
are defined in meaning and card column position in Table 1 and illustrated 
in Figures l-3. 


COMPARISON OF COMPUTER RESULTS WITH MANUAL CALCULATIONS 


Impulse - PSI - MS 


— 


Problem Manual __——_—s Computer Difference % Difference 
AS 1000 958 42 4.2 
A6 8500 8651 151 1.78 
i _A8 117 Liz 5 4.24 | 


As indicated in Tabie 2, the difference between the manual and 
computer calculations is only a few percent. This difference may be 
accounted for by the smoothing of the impulse curves contained in the 
manual. The computers results should »e more accurate than the 
manual charts since the calculation are direct and numerical inter- 
polation is used instead of visual. 


Appendix A shows the details of the manual and computer calculations. 
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These procedures should be followed in filling out the input form. 


1. Fill in Title or Problem identification Card. 
. 2. Sketch the structure and charge location as shown in Figure l 
and 2. 


e 3. Enter RA, W, H, L. h, 1. (Units are in feet and ibs; a decimal 
| point must be supplied after each number.) 


4. Enter a"l" or "zero" (no decimal point) in the appropriate 
column of the reflection code of the input sheet. If the re- 
' flection surface (floor, rocf, left side wall or right side 
; wall) is present enter a 1, otherwise indicate the absence of 
i the surface by a zero. Sketching the cubicle and unfolding the 
7 E view (Figure 3) will help determine the reflecting surface. 
‘k ' Detailed instructions follow. 
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; . Specifying Reflection Code t 


In perform..ng its calculations, the program considers the effcct of 
reflection of tne original blast impulse from surfaces at right angles and 
; in contact with the wall whose impulse load is being computed. Figure 3 
illustrates a method for determining the reflecting surfaces of a three- 
wall cubicle with roof. Figur2 3(a) shows how to specify the reflection 
/- code on the input sheet (Figure 4) when it is desired to calculat: the 
impulse load on the back wall. By unfolding the walls of the cubicle, ; 
keeping the back wall in the center, it is noted that there are four re- 
flecting surfaces at right angles and in contact with the back wall: 
the floor, roof and two side walls. There the presence of these reflecting 
surfaces are indicated on the input sheet by putting a Ll in each cf the 
<a columns labeled Floor, Roof, Left Side, Right Side. The reflection code 
of Figure 4 would be 1111. Sample Problcan A6 Q@ppendix A) has the same 


configuration of Figure 4 except that it lacks a roof; its reflection t 
code is 1011. 


Sometimes the wall in question is not the back wall. In order to 
use the input form, the cubicle should be rotated so that the wall 
/ occupies the same position as the back wall. The solution will be the 
" § same as a back wall problem. 
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In Sample Problem A8, (Appendix A) it is required to calculate the 
impulse load on the roof of a cubicle. The cubicle is rotated 90° so 
that the roof becomes the back wall, the floor the front wail, and the 
front wall the roof. This configuration will then be the same as 
Figure 3(a). The input sheet should then be filled out accordingly. 


Figure 3(b) shows an unfolded view of the same cudicle used as an 
aid in calculating the impulse on a side wall. The three reflecting 
surfaces are the roof, back wall and floor. In specifying the input code, 
the cubicle is rotated 90° so that the side wall is treated as a back 
wall and vice versa. The reflection code of Figure 4 (input sheet).would be 1110. 
Sample Problem A5 (Appendix A) shows the calculation of the {mpulse 
loading on a side wall of a cubicle slmilar to Figure 4 except that 
Sample Problem AS lacks a roof; its reflection code is 1010. 
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